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ABSTRACT 

A complete understanding of the origin of the prestellar core mass function (CMF) is crucial. Two major features of the prestellar CMF are: 1) 
a broad peak below 1 Mq, presumably corresponding to a mean gravitational fragmentation scale, and 2) a characteristic power-law slope, very 
similar to the Salpeter slope of the stellar initial mass function (IMF) at the high-mass end. While recent Herschel observations have shown that 
the peak of the prestellar CMF is close to the thermal Jeans mass in marginally supercritical filaments, the origin of the power-law tail of the 
CMF/JMF at the high-mass end is less clear. Inutsuka (2001) proposed a theoretical scenario in which the origin of the power-law tail can be 
understood as resulting from the growth of an initial spectrum of density perturbations seeded along the long axis of star-forming filaments by 
interstellar turbulence. Here, we report the statistical properties of the line-mass fluctuations of filaments in the Pipe, Taurus, and IC5146 molecular 
clouds observed with Herschel for a sample of subcritical or marginally supercritical filaments using a 1-D power spectrum analysis. The observed 
filament power spectra were fitted by a power-law function (F’tiuc('S) ■?“) after removing the effect of beam convolution at small scales. A 
Gaussian-like distribution of power-spectrum slopes was found, centered at dcorr = -1.6 ± 0.3. The characteristic index of the observed power 
spectra is close to that of the one-dimensional velocity power spectrum generated by subsonic Kolomogorov turbulence (-1.67). Given the errors, 
the measured power-spectrum slope is also marginally consistent with the power spectrum index of -2 expected for supersonic compressible 
turbulence. With such a power spectrum of initial line-mass fluctuations, Inutsuka’s model would yield a mass function of collapsed objects along 
filaments approaching dN/dM oc yi^-2.3±o.i high-mass end (very close to the Salpeter power law) after a few free-fall times. An empirical 

correlation, P^-^(sq) oc *, was also found between the amplitude of each filament power spectrum P{sf) and the mean column density 

along the filament (Ahj). Finally, the dispersion of line-mass fluctuations along each filament was found to scale with the physical length 
L of the filament, roughly as oc Lf-''. Overall, our results are consistent with the suggestion that the bulk of the CMF/IMF results from the 
gravitational fragmentation of filaments. 
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1. Introduction 

Understanding the origin of the stellar initial mass function 
(IM F) is a fundamenta l open problem in modern astrophysics 
(e.g. lOffner et al.ll2014l for a recent review). Since the end of 
the 1990s, several observational studies of prestellar dense cores 
in nearby molecular clouds have found a strong link between 
the prestellar core m a ss function (CMF) and the IMF (e.g., 
see | Motte et alJ Il998l : lAlves et ^ 120071 : iKonvves et alJ l20lH 
120 li T suggesting that the IMF is at least partly the result of 
the core formation process. Theories of the CMF/IMF based on 
gravo-turbulent fragmentation (e.g. iPadoan & Nordlundl l2002t 


iHennebelle & Chabrieti20(MlHoDkinsl2012 ) are consistent with 


Herschel observations fe.s.. lAndre et al.1 

201C 

Molinari et al. 

201C 

: lArzoumanian et alj|201 It iHill et alJ 

2011 

Konvves et al. 


this view but do not account for the fact that most cores/stars 
appear to form within interstellar filaments. Indeed, recent 


supercritical filaments with masses per unit length Mune ap¬ 
proaching the critical line mass of nearly isot hermal ~ 10 K 
gas cylinders MiinR.mt = 2c?/G ~ 16 Mo/pc (IOstrikeilll964 
llnutsuka & Mivamalll997l) may be responsible for the peak of 
the prestellar CMF at ~ 0.6 Mn, as observed in the Aquila clou d 
complex for example (lAndre et al.l 120141 : iKonvves et al.llT015l) . 
Indeed, this idea is consistent with the view that the peak of the 
IMF is related to the typical Jeans mass in star-forming clouds 
(lLarsonl|1985[) . 

In gravo-turbulent fragmentation theories of the CMF/IMF, 
the peak of the IMF results from a combination of thermal 
physics, setting the mean thermal Jeans mass in the parent 
clo ud, and turbulence effects, through the turbulence Mach num¬ 

ber (Hennebelle & Chabrierl2008l:lHoDkinsl2012l:ICha brier et aF 



201-51) emphasize the role of interstellar filaments in star-forming 
clouds and support a paradigm for star formation in which the 
formation of ~ 0.1 pc-wide filaments and the subsequent frag¬ 
mentation of the densest filaments into prestellar cores repre¬ 
sent t wo key steps in the star formation process (cf. lAndre et alJ 
|2014|) . The results of the Herschel Gould Belt survey (HGBS) 
further suggest that gravitational fragmentation of marginally 


In the filamentary picture proposed bv lAndre et al 

a characteristic thermal Jeans mass results from the ex¬ 
istence of a critical line mass for filaments (which depends only 
on gas temperature) and from the c haracteristic filament widt h 
~ 0.1 pc measured with Herschel (lArzoumanian et akl 1201 ll) . 
which is close to the sonic scale of turbulence in low-density 
molecular gas. Taken together, this sets a (column) density 
threshold for prestellar core and star forrnation as observed in 
nearby clouds with Herschel dAndre et al.ll2010t IKonvves et al.1 
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Fig. 1. Herschel 2 50 jim surface brightness image of the Pipe molecular cloud observed as a part of the HGBS key project (cf. 
iPeretto et alJl2Q12]l . The image is displayed a t the diffraction limited 250 jim beam resolution of 18"2. The 43 filaments traced with 
the DisPerSE algorithm (ISousbie et al.ll20111) are highlighted in white. A more detailed view of the filament enclosed by the black 
dashed rectangle is shown in Fig. 


1201 5h and Spitzer (iHeiderman et al.ir2010t iLada et aTll2010l) . re¬ 
spectively. 

Although the thermal Jeans mass within marginally super¬ 
critical filaments may account for the peak of the prestellar 
CMF, pure thermal fragmentation of filaments at the thresh¬ 
old of gravitational instability is expected to lead to a sharply 
peaked CMF (such as a delta function) and can hardly explain 
the Salpeter pow er-law regime of the IMF at the high-mass end. 
Ilnutsukal(l200lh nevertheless suggested that a Salpeter-like CMF 
can quickly develop within filaments provided that turbulence 
has generated the appropriate power spectrum of initial density 
fluctuations in the first place. More specifically. Ilnutsuka d200 ih 
found that the hierarchical structural property of the perturbed 
density field along filaments may statistically produce a popu¬ 
lation of cores with a power-law mass function. Inutsuka’s ap¬ 
proach consisted of counting the distribution of “isolated” col¬ 
lapsed regions of mass scale M above a critical density using the 
Press -Schechter formalism dPress & SchechteJl 19741: IJedamzil3 
Il995h . where the statistics of such regions depend on the log¬ 
arithmic slope of the power spectrum characterizing the initial 
fluctuating density field. Inutsuka (2001) showed that the re¬ 
sulting core mass distribution converges toward a Salpeter-like 
mass function dNIdM oc when the power spectrum of the 
field of initial line-mass fluctuations approaches P{s) oc 
The latter power spectrum was adopted on a purely ad-hoc ba¬ 
sis, however, and the true statistical nature of the density fluc¬ 
tuations along filaments remained to be quantified observation- 
ally. Another important assumption of Inutsuka’s model was that 
the initial line-mass fluctuations along filaments could be repre¬ 
sented by a Gaussian random distribution. 

In this paper, we exploit the unprecedented sensitivity and 
resolution of Herschel submillimeter continuum images as well 
as the proximity of the molecular clouds targeted as part of the 
HGBS key project to characterize, for the first time, the sta¬ 
tistical properties of the line-mass fluctuations along interstel¬ 
lar filaments by analyzing their 1 -D power spectrum. Our main 
goal is to constrain observationally the validity of the theoretical 


model outlined by llnutsukal (1200ll) for the origin of a Salpeter- 
like prestellar CMF above ~ 1 Mq. 


2. Herschel observations and column density maps 


The three target fields that we used in the present analysis were 
all ob served with the Hersche^ space ob servatory (Pilbratt et al.l 
1201 oh as part of the HGBS key project (I Andre et ^ 201 Oh . and 


cover surface areas of ~ 1?5 x 1?5, ~ 6° x 2?5, and ~ 1?6 x F6, 
in the ICS 146, Pipe, and Taurus molecular clouds, respectively 
(lArzoumanian et ^ 1201 ll : iPeretto et al.ll2012t iPalmeirim et al.l 
l2013h . These three regions were mapped at a scanning speed 
of 60"s ' in parallel mode si multaneously at fiv e Herschel 
wavelengths using th e SPIRE (iGriffin et al J 1201 Oh and PACS 
(IPoglitsch et al.ll201(^ photometric cameras. The data were re¬ 
duced using HIPE version 7.0. For the SPIRE data reduc¬ 
tion, we used modified pipeline scripts. Observations during 
the turnaround of the telescope were included, and a destriper 
module with a zero-order polynomial baseline was applied. The 
default ‘naive’ mapper was used to produce the final maps. 
For the PACS data, we applied the standard HIPE data reduc¬ 
tion pipeline up to level 1, with improved calibration. Further 
processing of the data, such as subtraction of (thermal and 
non-thermal) low-frequency noise a nd map project ion was per¬ 
formed with Scanamorphos vll (iRousse J l2013h . Note that 
the Scannamorphos map-maker avoids any high-pass filtering, 
which is crucial for preserving extended emission. We adjusted 
the zero-point values of individual Herschel images based on 
cross-cor relations of the Her schel data with IRAS and Planck 
data (cf. [Bernard et akl l2010l) . Figure [1] shows an example of 
250 pm surface brightness image of the Pipe Nebula at the native 
resolution of 1872. 

For the present analysis, we needed high-resolution col¬ 
umn density maps for recovering small-scale fluctuations which 


* Herschel is an ESA space observatory with science instruments pro¬ 
vided by European-led Principal Investigator consortia and with impor¬ 
tant participation from NASA. 
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would otherwise be smeared due to beam convolution. We em¬ 
ployed two methods for generating such column density maps. 
First, we generated column density maps at an effective reso¬ 
lution corresponding to the SPIRE 250 rtm beam fol l owing the 
method described in Appendix A of iPalmeirim et al.l (l2013h . In 
this method, large-scale information is obtained from a ‘stan¬ 
dard’ Herschel column density map constructed at the SPIRE 
500 fim resolution, while hner details are recovered by con¬ 
structing more approximate column density maps at the reso¬ 
lution of SPIRE 350 jum and SPIRE 250 /rm data, respectively. 
Eollowing the steps described in Palmeirim et al., we fitted a 
single-temperature modihed blackbody to the observed data on 
a pixel by pixel basis, adopting a wavelength-dependent dust 
opacity of the form k,i - 0.1 (/i/300 jum)’^ cm^ per g (of gas 
+ dust) with an emissivity index — 2 (cf. [Hildebrandlll983t 
IRov et al]l2014l) . In the second, alternative method, we directly 
converted the surface brightness maps observed at 250 fim I 250 
into approximate column density maps by using the relation 
A^H2='l25o/(B25o[7’d,sED]'f250/^H2'«H), where/ 2 h 2 = 2.8 is the mean 
molecular weight and Fd.sEO is the SED dust temperature ob¬ 
tained by htting a modified blackbody to the SED observed be¬ 
tween 160 ijm and 500 fim with Herschel data toward each line 
of sight. Eor the Pipe Nebula filament shown in Pig. |2] for in¬ 
stance, the median dust temperature along the filament crest is 
Fd.SED = 13.5 K. The advantage of the first method is that it 
produces a more accurate column density map at the 18.2" res¬ 
olution of the SPIRE 250 data. The disadvantage, however, 
is that the resulting map is signihcantly noisier than the SPIRE 
250 fim map and the point spread function (PSP) is more difficult 
to characterize due to the various steps involved in the process¬ 
ing. In contrast, the second method which simply produces mod¬ 
ified 250 yum maps, the effective PSP is exactly the same as that 
of the SPIRE 250 fim photometer. Since an accurate knowledge 
of the beam is important in the present power-spectrum study, we 
report results based on the latter method in the main body of the 
paper. The results of a power-spectrum analysis performed using 
high-resolution column density maps produced with the former 
method are given in Appendix B. They conhrm the robustness 
of our conclusions. 


3. Filament identification 


Eor the purpose of detecting coherently elongated filamentary 
structures in the column density maps as well as the modified 
250 fim images, we employed the DisPerSE algorithm which 
traces the crests of (segments of) “topological” filaments by con¬ 
necting saddl^_goii^tomaxima fol lowing the gradient in an 
image ([Sousbie et al.l 1201 ll see also lArzoumanian et al. IUMI 
iHill et ai.ll201 ll for details on practical applications of DisPerSE 
on Herschel images). It is important to point out here that 
DisPerSE was used solely for tracing the filament crests, while 
the hlament longitudinal profiles were measured directly on the 
original images. Our filament sample consists o f 28 hlamen ts in 
the Taurus cloud (adopted distance d - 140 r)c: lElias Ill978h. 42 
filame nts in the Pipe nebula region (d = 145 pc; Alves et ^ 
2007h. and 36 f ilaments in the IC5146 cloud {d - 460 pc; 


Lada et aDll999h . Combining the three regions, we identihed 
a total of 106 filaments of various lengths above a 5 cr persis¬ 
tence level of ~ 10^* cm“^. For a proper power-spectrum analy¬ 
sis which was our ultimate goal, however, we had to select mod¬ 
erately long filaments to have an adequate number of Fourier 
modes at small spatial frequencies (see Sect.|5]below). Our final 
analysis was thus based on the subset of 80 filaments which are 
longer than 5'5 (~ 18 beams). 
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Fig. 2. Close-up view of the sub-region marked by a black 
dashed rectangle in Fig. [1] showing (a) an example of a nearly 
critical filament in the Pipe cloud at an adopted distance d - 
145 pc, and (b) the line-mass/column density fluctuations along 
the long axis of the same filament, both at the 18'.'2 (HPBW) 
resolution of the SPIRE 250 data. The black curve shown 
in (a) marks t he crest of the filam ent as traced by the DisPerSE 
algorithm (cf. ISousbie et al.|[201 ll) . In b) the right ordinate axis 
shows the net H 2 -column density fluctuations along the filament. 
Equivalent line-mass values (assuming a characteristic filament 
width of 0.1 pc) are indicated on the left ordinate axis. A constant 
background column density level was subtracted before deriving 
the line mass fluctuations. The dotted horizontal line marks the 
mean line mass of the filament, (Miine) ~ 15 Mo/pc, which is 
only slightly lower the critical line mass of 16.3 Mo/pc for an 
isothermal cylinder at 10 K. The median relative amplitude of 
the line-mass fluctuations is ~ 0.07, within the linear 

(j^line) 

perturbation regime. The maximum relative amplitude is ~ 0.3 
« 1 . 


4. Column density/line mass fluctuations along 
selected filaments 

Eigure |2| shows an example of a longitudinal profile of column 
density or line-mass fluctuations along a marginally critical fil¬ 
ament in the Pipe molecular cloud. As the mass per unit length 
or l ine mass is a f undamental variable in cylindrical geometry 
(cf. Ilnutsukall200ll) . we converted the observed column density 
fluctuations to line-mass fluctuations by mak ing use of the quasi- 
univer sal filament width ~ 0.1 pc found bv lArzoumanian et al.l 
(1201 ih . The line-mass fluctuations corresponding to the exam¬ 
ple filament displayed in Pig. are shown in Fig. |2]5, while 
the equivalent column density fluctuations are quantified on the 
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Fig. 3. Normalized histogram showing the probability distribu¬ 
tion of relative line-mass fluctuations along the long axes of the 
subset of 67 filaments which are entirely subcritical. The red 
curve shows the Gaussian best fit to the observed distribution, 
which has a mean of ~ 0 and a standard deviation of ~ 0.09. 


right ordinate axis of Fig. |2jr. To estimate fluctuations intrinsic 
to the selected filaments themselves, a mean line-of-sight back¬ 
ground column density was subtracted from the observed data 
for each held. The dotted horizontal line in Fig. |2j) shows the 
mean line-mass of the filament (~ 15.5 Mo/pc), allowing one to 
evaluate the relative strength of the line-mass fluctuations about 
the mean. In this example, as well as for all other subcritical 
filaments in our sample, the line-mass fluctuations are small in 
relative amplitude, i.e., < 0 . 1 , implying that the per¬ 

turbation modes are in the linear regime (see also Sect.j^below). 

Figure|3]shows the distribution of normalized line-mass fluc¬ 
tuations, (5 = ^ filaments of our sample 

which are subcritical over their entire length. It can be seen in 
Fig.[2that the distribution of line-mass fluctuations for subcriti¬ 
cal filaments is nearly Gaussian, which is an important assump¬ 
tion of the CMF/IMF model calculations presented by llnutsukal 
(1200 ih . It can also be seen that the density/line-mass fluctuations 
along subcritical filaments are typically less than 10 %, i.e., are 
in the linear regime. 

The power spectrum of the fluctuations observed along the 
filament displayed in Fig. is shown in Fig. |4] and discussed in 
Sect. |5]below. Panels (a) and (b) of Fig. |5] and Fig. | 6 ] are similar 
to Fig. |2]but show two extreme examples of a subcritical fila¬ 
ment (with < Miine > ~ 3 Mo/pc << Miine,crit) and a supercritical 
filament (with < M^ne > ~ 30 Mo/pc > Miine,cnt) in the IC5146 
and Taurus clouds, respectively. 

Note that the column density/line mass longitudinal profiles 
used in this paper (see, e.g., Figs.|2}r &|23 &| 6 ]r) have an effective 
18'.'2-beam resolution. The pixel resolution of the corresponding 
data is 6'.'0, close to the Nyquist sampling limit {AjlD ~7'.'4) 
of SPIRE 250-/zm observations where D - 3.5m is the diame¬ 
ter of the primary mirror of the Herschel telescope. When de¬ 
riving physical properties from the power spectra of individual 
filaments, we took special care to exclude angular frequencies 
higher than the Nyquist critical frequency DjA ~ 1/(2 x 7.4") ~ 
0.068 arcsec^' ~ 4.0arcmin“'. 
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Fig. 4. Power spectrum of the line-mass fluctuations observed 
along the Pipe filament shown in Fig. as a function of spa¬ 
tial frequency s (bottom x-axis) or angular frequency s (top 
x-axis). In both panels, the black plus symbols show the ob¬ 
served power spectrum, PobsC*)- In panel a), the red dots show 
the power spectrum, PohAs) - Pn, obtained after subtracting a 
white noise power spectrum level p^timated 1 x 10 “^ Mo^/pc, 
marked by the horizontal red line and estimated from the me¬ 
dian value of Pobs('S) in the 3.9^.2 arcmin^’ angular frequency 
range; the blue triangles show a similar power spectrum after 
subtracting the instrument noise power spectrum level ~ 

4.0x10“^ M 0 ^/pc, marked by the horizontal blue line and cor¬ 
responding to the instrument noise level Icr ~ 1 MJy/sr in our 
SPIRE 250 ;um maps according to HSpot^. In panel b), the cyan 
dots show the noise-subtracted and beam-corrected power spec¬ 
trum, Ptrne('S) = (PobsC*) - PN/Zybeam- The Vertical dotted line in 
both panels marks the FWHM of the beam power spectrum at 
250 yum (S Si 2 arcmin^' - see Eig. lA.ll l. which is also the high¬ 
est frequency data point used to fit a power-law function. The 
vertical dot-dashed line is the Nyquist angular frequency (/1/2D) 
for SPIRE 250 jjm data. The power-law fits to the power spectra 
Pobs(s) and PtrueC*) have logarithmic slopes cTobs = - 2.1 + 0.2 and 
Qftme = -1.6 + 0.2, respectively. (Considering only angular fre¬ 
quencies up to S = 1.5 arcmin^*, the best power-law fit to Ptrue('S) 
has a slope atme = -1.7 ± 0.3.) 

5. Power spectrum analysis of the filaments 

The normalized power spectrum P{s) of a ID spatially-varying 
held 7(0 is proportional to the square of the Eourier amplitude 
of the signal, and mathematically the relationship in the one¬ 
dimensional case is: 

P{S) = i|7(^)|2, (1) 

where s denotes spatial frequency (and s denotes angular fre¬ 
quency), I(s) - J /(/) e'ZOTi/ j-jjg Fourier transform of 1(1), 

and L - fdl is the total length over which the signal is mea¬ 
sured. In the present case, it is convenient to take /(I) to be the 
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Fig. 5. Panels (a) and (b) are similar to Fig. |2^ and Fig. |2]r but 
for an example of a thermally subcritical filament in the IC5146 
molecular cloud (d~ 460 pc). The mean line mass of the fila¬ 
ment is < Miine > ~ 3 Mo/pc << Miine.crit. as shown by the 
dotted horizontal line in panel (b). The median relative ampli¬ 
tude of the line-mass fluctuations is ~ 0.06, within 

the linear perturbation regime. Panel (c) is similar to Fig. Ubut 
for the IC5146 filament. The power-law fits to the power spec¬ 
tra Pohs(s) and Pti-ue(‘^) have logarithmic slopes -1.8 + 0.2 and 
-1.2 + 0.4, respectively. (Considering only angular frequencies 
up to 5 = 1.5 arcmin" \ the best power-law fit to Ptrue(‘^) has a 
slope ati-ue = -1.1 ± 0.3.) 
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Fig. 6. Similar to Fig. |5] but for the thermally supercritical 
filament B211/B 213 in the Taurus cloud at d~ 140 pc (cf. 
iPalmeirim et alJl2013f) . The mean line mass of the filament is 
< Miine >~ 30 Mo/pc >> Miine,crit, as shown by the dotted hori¬ 
zontal line in panel (b). (Here, < Mune > was obtained by mul¬ 
tiplying the mean background-subt racted column density by a 
characteristic 0.1 pc width, whereas [Palmeirim et al.l (1201 3li ob¬ 
tained a value of 54 Mo/pc by integration over the full radial 
column density profile of the filament.) The median relative am¬ 
plitude of the line-mass fluctuations is ~ 0.4. Note 

the spike features corresponding to the positions of dense cores 
along the filament and capturing non-linear density perturba¬ 
tions. In panel (c), the power-law fits to the power spectra Pobsi^) 
and Ptrue('S) have logarithmic slopes -2.2 + 0.2 and -1.9 + 0.3, 
respectively. (The best power-law fit to Pti-ue('S) is unchanged if 
only angular frequencies up to s = 1.5 arcmin ' are considered.) 


held of line-mass fluctuations (in Mo/pc) along a given filament 
and to express the offset along the filament crest, I, in units of pc. 
Then, the normalized power spectrum P(s) has units of My-pc. 

The observed signal after convolution with the telescope 
beam, B, is given by /obs(0 = ftrue(0 * B + N{J), where N(l) 
is a noise term arising from small-scale fluctuations primar¬ 
ily due to instrument noise. The Fourier transform of /obs(0 is 
fobs('S) = hmeis) X B{s) + N{s). Using Eq. ([TJ and the foregoing 
relation, it can be seen that the total power spectrum observed 
along the filament axis is composed of the true power spectrum, 
Ptrue(s), due to liue-mass variations along the filament modified 
by the power spectrum of the beam, ybeam('S). plus a white noise 


power spectrurrQ term, Pbi(s), and can be expressed as: 

Fobs(‘^) — Ft];i^[e(5) X ybeam(‘^) Fbj(.s). (2) 


^ The SPIRE bolometers have a measured temporal stability S 250 s, 
which for Herschel parallel-mode observations taken at a scanning 
speed of 60" s“’ translates to a 1// noise ‘knee’ at spatial frequen- 
cies corresponding to angular scales larger than ~ 4° (see Fig. 2 of 
iPascale et alj|201]]) . The angular scale of the longest filament consid¬ 
ered here ('- 2°) is significantly smaller than both the angular scale of 
the 1 // knee and the maximum scan length of the corresponding scan- 
map observations. We can therefore safely assume that the instrumental 
noise behaves like a white noise over the entire range of spatial frequen¬ 
cies considered in the paper. 
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Fig. 7. Distributions of power spectrum slopes measured before 
beam correction (dashed histogram) and after beam correction 
(solid histogram) for our selected sample of 80 filaments. Best- 
fit Gaussian curves to the two observed distributions are over¬ 
plotted. The two distributions are centered on dobs = -2.0 + 0.2 
and dcoiT = -1.6 + 0.3 for the uncorrected and beam-corrected 
power spectra, respectively. 

In deriving Eq. (2), we assumed that the noise and the filament 
line-mass fluctuations are completely uncorrelated, resulting in a 
vanishing cross-power spectrum term. To estimate the amplitude 
of P^is), we computed the median of the power spectrum in a 
narrow band of angular frequencies, 3.9 arcmin ’ < s < 4.2 
arcmin *, centered about the Nyquist angular frequency of 4.0 
arcmin"*. 

In practice, to recover PtmeC*). we first computed Pobs)*) 
directly from the data, then subtracted a constant noise power 
level, and Anally corrected the result for the beam convolution 
effect by dividing by 7beam('S)- Within the relevant range of an¬ 
gular frequencies (5 < FWHM of the beam power spectrum 
~ 2 arcmin"*), both Pobsi^) and Ptrne('S) can be approximated and 
fitted by a power-law function: 

Pobs/tme(i) ~ P(io)(i/io)“*'''"% (3) 

where P{so) is the normalizing amplitude of the power spectrum 
at spatial frequency sq and a is the power-law index of the fitted 
power spectrum. 

The results of such a power spectrum analysis are shown 
in Fig. |4] for the marginally critical filament of Fig. |3 The 
plus symbols in Fig. |4] correspond to the observed power spec¬ 
trum of the example filament, Pobs('S). prior to noise subtrac¬ 
tion and beam correction. The upper panel of Fig. |4] explic¬ 
itly shows the effect of noise level subtraction. The horizon¬ 
tal red line marks the estimated power level due to instrument 
noise, p^stimated _ jxlO"'* Mg^/pc, derived following the above 
recipe. For comparison, the blue horizontal line shows the noise 
power spectrum level (P^^*’°') corresponding to the rms instru¬ 
ment noise level (i.e., Icr ~ 1 MJy/sr or ~ 10 mJy/beam) ex¬ 
pected for SPIRE 250 pm data corresponding to two orthog¬ 
onal scans taken at 60"s"* in parallel-mode observations ac¬ 
cording to the Herschel Observation Planning Tool (HSpotfl 
The noise power spectrum level was determined using 

Parseval’s theorem which relates the variance of the instrumen¬ 
tal noise to the integral of its power spectrum: CTMiineHSpot “ 
pmax _ 2 ^ 7’^^*’°', where the integral is taken over 

the relevant spatial frequency range given the maximum angular 
frequency Sniax = 4.0 arcmin * sampled by our SPIRE 250 pm 
data. For a 0.1-pc wide filament with a median dust temperature 

^ http:// herschel.esac.esa.int/Docs/HSPOT/html/hspot-help.html 


of 13.5 K as measured along the filament crest in Fig. |2] a Icr 
rms level of ~ 1 MJy/sr at 250 pm corresponds to a standard 
deviation CTMime.HSpot ~ 0.09/V/o/pc in line-mass units. Given 
the distance d - 145 pc to the filament of Fig. |2 this yields 
= *^Mme,HSpot/(2 imax) ~ 4.0 X 10"^/w|/pc in the Same 
unit as that used on the y-axis of Fig. |4] 

The red dots and overplotted blue triangles in Fig. |4^ show 
the power spectrum data points after subtracting the estimated 
noise spectrum level and the expected instrument noise 

spectrum level respectively. It is apparent from Fig. |4^ 

that the departure of the noise-subtracted data points from the 
observed filament power spectrum is negligible in the angu¬ 
lar frequency range s < 2 arcmin"* where we fit a power- 
law slope. Therefore, even though our estimated noise spectrum 
level, p^timated^ likely overestimates the actual noise spectrum 
level, the step of noise subtraction has effectively no influence 
on the estimated power-law slope of the observed power spec¬ 
trum. 

The beam-corrected power spectrum, Ptmei^), shown by the 
cyan solid dots in Fig. I4J) was obtained by dividing the noise- 
subtracted power spectrum by the 250 pm SPIRE beam power 
spectrum (ybeam)- The derivation of the SPIRE beam power spec¬ 
trum is described in AnpendixlAl The Pobs)*) spectrum has more 
power at low spatial frequencies s (large scales) and decreases 
roughly as a power law toward high spatial frequencies until it 
merges with the noise power spectum level at a spatial frequency 
close to s ~ 100 pc"* (i.e., angular frequency ~ 4.0 arcmin"*). 
The beam-corrected power spectrum, Ptmeis), has a slightly shal¬ 
lower slope than Pohsis)- At spatial frequencies higher than the 
FWHM of the beam power spectrum shown by the dotted ver¬ 
tical line at i ~ 50 pc"* (angular frequency s ~ 2 arcmin"*) in 
Fig. a Ptrue('S) is overcorrected due to amplification of residual 
noise when dividing by the beam power spectrum. This unphys¬ 
ical rise of Pti-ue)*) (see Fig. |4j)) at very small angular frequen¬ 
cies is a generic feature of beam-corrected power spectra and 
is w ell documented in earlier Herschel studies of cirr us noise 
te.g. iMartin et al.ll2010( iMiville-Deschenes et alJl2010ll . To get 
around this problem, it is customary to avoid angular frequen¬ 
cies greater than the FWHM of the beam power spectrum, i .e., 
5 > 2 arcmin * at 250 pm (cf. IMiville-Deschenes et al1l2010l) . a 
compromise that we also adopt in this paper. Upon visually in¬ 
specting Fig.|4j3, there may still be an indication of a slight over¬ 
correction of Ptmsis) in the frequency range 1.5 arcmin"* < S < 2 
arcmin"*. Therefore, we also report the results of more conserva¬ 
tive power-law fits obtained by fitting PtraeC*) only up to s = 1.5 
arcmin"*. Finally, it is common practice to avoid very small spa¬ 
tial frequencies when fitting power spectrum slopes in the case 
of 2D images, due to reasons such as edge effects or large scan¬ 
ning lengths. In principle, in the case of a filament much smaller 
than the size of the mapped region, there should be no techni¬ 
cal difficulty in reconstructing features at all angular scales with 
Fourier modes. To be on the safe side, we nevertheless excluded 
the smallest non-zero angular frequency data point when fitting 
a power law to the power spectrum of each filamenfl- 

For the Pipe filament shown in Fig. |4l the beam-corrected 
power spectrum has a best-fit slope atrue=-l-6 + 0.2 and an 
amplitude P(so = 10 pc"*)= 9.5x10"^ /Wo^/pc. For compari¬ 
son, the observed power spectrum has a somewhat steeper best- 
fit slope aobs=-2.1 + 0.2 and an amplitude P(so - 10 pc"*)= 
6x10"^ /Wo^/pc. Two other examples of filament power spectra 

Note that including the lowest s data point when fitting the power 
spectrum slope would change the overall distribution of slopes in Fig.|7] 
by less than one standard deviation. 
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are shown in Fig.|5}; and Fig. |6j; for the ICS 146 subcritical fila¬ 
ment of Fig. |5^ and the Taurus supercritical filament of Fig. |6^, 
respectively. 


N„2 [10^' cm“^] 



<M|;ne> [^o/P*^] 

Fig. 8. Square root of filament power spectrum amplitude at so 
= 10 pc ' versus mean line mass of each filament. The solid 
line shows the best-fit power-law The cor¬ 

responding mean Ahj- column density level is shown by the up¬ 
per coordinate axis. The vertical dashed line separates subcrit¬ 
ical filaments (cyan filled circles on the left) from supercritical 
hlaments (blue filled circles on the right). 


5.1. A characteristic power spectrum siope for fiiaments 

To reach statistically significant conclusions on the typical 
power spectrum slope, we considered a sample of 80 subcritical 
(or marginally supercritical) filaments belonging to three distinct 
nearby molecular clouds (see Sect.O. Since the reliability of the 
derived power spectrum properties also depends on the dynamic 
range covered in spatial frequency space, we excluded from our 
present analysis filaments which were smaller than 5f5 in the 
plane of sky. The resulting distribution of power spectrum slopes 
is shown in Fig.|2] The two histograms shown by the dashed and 
solid lines represent the distributions obtained for the observed 
[F’obs(i)] and beam-corrected [F’tnje(i)] filament power spectra, 
respectively. Both histograms are well fitted by Gaussian distri¬ 
butions (see Fig.|7]i. 

The mean power-law index dobs = -2.0 + 0.2 measured for 
the observed power spectra is slightly steeper than the mean 
power-law index dcon = -1.6 + 0.3 obtained at s < 2 arcmin”' 
after correcting for the beam convolution effect. Fitting Paueis) 
with a power law over the more conservative range of angular 
frequencies i < 1.5 arcmin ’ yields essentially the same index 
within errors, dcorr = -1.8 + 0.5. The distribution of beam- 
corrected power spectrum slopes has a larger dispersion than the 
distribution of uncorrected slopes, due to the propagation of un¬ 
certainties resulting from the amplification of residual noise at 
small scales. 

5.2. Correiation between P{so) and mean coiumn density 

Figure 0 shows a positive correlation between the amplitude of 
the line-mass fluctuations, expressed by P'^^(so), and the mean 
column density of the filaments. In other words, filaments with 
higher mean column densities tend to have higher power spec¬ 
trum amplitudes, i.e., stronger column density fluctuations. The 
red straight line in Fig.[8]shows the best-fit power law to the data 
points, which has an index of 1.4 + 0.1. Physically, this empir¬ 
ical relation means that the rms fluctuations of the perturbative 



Length [pc] 


Fig. 9. Dispersion of line-mass fluctuations ctmh,,^ versus pro¬ 
jected length L for the filaments in our sample. The solid line 
shows the best-fit power-law ctmu,,, ^ to the data points. For 
comparison, the dashed line shows a power-law scaling equiv¬ 
alent to the well-kno wn linewidth- size relation in molecular 
clouds cTy oc (cf. lLarsonlFl981h . with an arbitary normal¬ 
ization. 


modes inside a filament increase with the mean line mass of the 
filament. A similar empirical correlation, with a somewhat dif¬ 
ferent power-law index (Pcimis {S v)^) exists between the 2D 
power spectrum amplitude of far-infrared cirrus (Pcimis) and the 
average surface brig htness {Sy) in far-infrared images of the sky 
(iGautier et al.lll99^ . 

5.3. Line-mass dispersion versus fiiament iength 

Figure |9] shows that the dispersion of line-mass fluctuations 
along the long axis of each filament, crM„„, is correlated with 
the physical length L of the filament. The best-fit power-law re¬ 
lation between ctmu,,, and L is crM„„ which is reminis¬ 

cent of the well-known power-law scaling between internal ve¬ 
locity dispersion an d region size in molecular clouds (cr^ oc ^) 
originally found bv iLarsonI (Il98lh . Such a similarity is not sur¬ 
prising since, in the subsonic turbulence regime which approx- 
imately holds for subcritic al or marginally subcritical filaments 
(lArzoumanian et al]l2013h . density fluctuations are expected to 
be directly proportional to velocity fluctuations. 


6. Concluding remarks 

Our one-dimensional power spectrum analysis along the axis 
of 80 nearby Herschel filaments shows that the longitudinal 
line-mass fluctuations along these filaments have a character¬ 
istic ID power spectrum slope dcon = -1.6 + 0.3. This re¬ 
sult is the first observational confirmation that the density/line- 
mass fluctuations along interstellar filaments have a characteris¬ 
tic power spectrum slope. The slope we find is remarkably close 
to the fiducial value atheoiy = -1 -5 adopted by llnutsukal (1200 ih 
in his attempt to explain the origin of Salpeter-like power-law 
tail of the CMF/IMF at the high-mass end {dNIdM oc 
Repeating Inutsuka’s analysis with the observationally-derived 
power spectrum slope (icon = -1 -6 found in the present stud}0 


^ Here, we adopt the mean, beam-corrected value dcorr =1.6 derived 
in Sect. 15. li as our best e s timate of the power spectrum slope. Using the 
prescription of llnutsuO ll200lh with the mean power-law index dobs = 
-2.0 measured before beam correction (see Fig.|7} instead would lead 
to a CMF approaching dN/dM oc after a few free-fall times. 
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would lead to a CMF approaching dNIdM oc at the high- 

mass end after a few free-fall times, which is even closer to a 
Salpeter mass function (dNIdM oc 

Interestingly, the measured power spectrum of density fluc¬ 
tuations along interstellar filaments shows a striking similar¬ 
ity to the energy spectrum or velocity power spectrum of in¬ 
compressible h ydrodynamic al turb ulence in ID, i.e., E{s) oc 
oc (iKolmo goroyll 194 flfl Note that, as long as self- 

grayity does not dominate, the continuity equation implies a di¬ 
rect proportionality between the Fourier modes of small density 
fluctuations and the modes of subsonic yelocity perturbations. 
Our measured power spectrum slope for the density fluctuations 
(iicorr = -1.6 + 0.3) is thus consistent with the fact that subcriti- 
cal (and marginally supercritical) filaments, such as the filaments 
considered in the present study, are obseryed to haye subsonic or 
at most transonic turbulent yeloc ity dispersions dHacar & Tafallal 
1201 It [Arzoumanian et al.l[2013h . 

The origin of the characteristic power spectrum slope of fil¬ 
aments may haye a close connection with the process generat¬ 
ing ’cirrus’ fluctuations of far-infrared surface brightness in the 
interstellar medium. The statistical an alysis of Gala c tic ‘ci rrus’ 
fluctuations in IRAS (2D) images by iGautier et al.l (Il992h has 
suggested a plausible role of interstellar turbulence in shaping 
these fluctuatio ns. The power sp ectrum of ‘cirrus’ fluctuations 
obtained from iRoy et al.l (1201 Ot) and recent Her schel images 
(iMartin et al.ll20icit Miyille-Deschenes et al.ll201^ has a power- 
law slope of -2.7, yery similar to the slope of -8/3 =» -2.67 
expected for the yelocity power spectrum of Kolmogoroy turbu¬ 
lence in 2D. The latter 2D slope translates to a power spectrum 
slope of -5/3 -1.67 in IDQ 

Our Herschel results support the following speculatiye pic¬ 
ture. Filamentary structures form in the cold interstellar medium 
as a result of a combination of compre ssion and shear in 
large-scale magneto - hydrodynamic flow s (IPadoan et a 


20011 


iBalsara et al.l 1200 It iHennebell^ l2013l llnutsuka et al.l l2015h . 
Supersonic turbulent motions dissipate quickly within dense fil¬ 
aments, forming more stable and coherent structures, and leay- 
ing an imprint of subsonic fluctuations along the filaments. 
Molecular line obseryations of the inner parts of subcritical or 
marginally supercriti cal filaments do show subsonic or transonic 
yeloc ity dispersions (iHacar & Tafallall201 ll [Arzoumanian et aP 
1201 3h . supporting the yiew that the dissipation of turbulence may 
play a role in the formation and eyolution of these structures. 
The density fluctuations along supercritical filaments are Jeans- 
unstable and grow by local grayitational instability, generating a 
distribution of prestellar cores. 

The obseryational yerification of the existence of a character¬ 
istic power spectrum slope for line-mass fluctuations along inter¬ 
stellar filaments has therefore opened an interesting ayenue for 
understanding an important aspect of the origin of the prestellar 
CMF. Our findings suggest that an initial spectrum of density 
perturbations seeded by interstellar turbulence may be a pre¬ 
requisite for generating a Salpeter-like CMF at the high-mass 
end. The present study combined with preyious Herschel results 
therefore proyides insight into the possible origin of two dis¬ 
tinct features of the CMF: the thermal grayitational fragmenta- 


® Given the uncertainties associated with beam correction (see 
Sect. |5j, the mean beam-corrected index dcorr = -1.6 ± 0.3 

is also marginally consistent with the energy spectrum E(s) oc 
tv expected for supersonic compressible turbulence (cf. 
iPadoan & Nordlundll2002l . and references therein). 

^ The general relation between the Kolmogorov energy spectrum 
E{k) oc Ic^!^ and the velocity power spectrum in n-D space is E(k) oc 
Pv(k). 


tion scale of filaments ma y set the peak of th e lognormal base 
of the prestellar CMF (e.g. lAndre et al.ll2()T^ . while the eyolu- 
tionary characteristics of the underlying line-mass perturbations 
may help t o set the power -law slope of the CMF/IMF at the high 
mass end (Ilnutsukall2001l and this paper). A final cayeat should 
be mentioned, howeyer. There is mounting eyidence that true 
high -mass analogs to low-mass prestellar cores may not exist 
(e.g. lMotte et alJ2007h and that massiye protostars may acquire 
the bulk of their final mass from much la rger scales than a sin¬ 
gle prestellar core (e.g. lPeretto et al.f2013h . Further work will be 
needed to eyaluate the relatiye importance of initial density fluc¬ 
tuations and l^e-scale accretion in generating the high-mass 
end ofthelMiff 
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Appendix A: Deriving the power spectrum of the 
SPiRE 250 jum beam 

The power spectrum of the SPIRE 250 /tm beam was derived 
from an empirical PSF image obtained by the SPIRE ICC from 
scan map dat^B of Neptune in four directions. There are two sets 
of empirical PSFs available, one gridded at the nominal pixel 
size for each SPIRE passband and the other at a higher pixel 
resolution of Cf-'h. We used the latter product, but the angular 
extent of the footprint was only about 6^0, corresponding to an 
angular frequency scale of ~ 0.3 arcmin“\ insufficient to correct 
for the beam convolution effect on large scales (i.e., small s). 
Therefore, to extend the angular frequency range of the beam 
power spectrum to lower s values, we embedded the beam map 
of Neptune inside a larger map 25f0 on a side, and then padded 
the pixels outside the central Neptune insert using the best-fit 2D 
Gaussian approximation to the beam. 

The power spectrum of the SPIRE 250 /im beam as a func¬ 
tion of angular frequency, ybeamCS), is shown in Eig. lA.ll The 
power spectrum asymptotically converges to unity at small an¬ 
gular frequencies but rapidly declines toward high angular fre¬ 
quencies (small angular scales). We caution that the SPIRE beam 
spectrum has a slight suppression of power around 0.2 arcmin^^ 
< S < 1.5 ar cmin~^ compared to a pu re Gaussian shape. This 
proble m led iMartin et alJ ( 1201 Oh and iMiville-Deschenes et alJ 
( 1201 Oh to approximate the beam power spectrum using higher- 
order polynomial corrections to the Gaussian fit. The effective 
full width at half maximum (FWHM) of the SPIRE 250 fim 
beam power spectrurrPl is about 2 arcmin^*, as shown by the 
dashed curve in Eig. IA.il 

Appendix B: Results from high-resolution column 
density maps 

Figure IB. II shows the power spectrum of the same filament al¬ 
ready displayed in Fig.|2]but derived from a high-resolution col- 

® The reduction of the beam map from 
scan-map data of Neptune is described at 
http://herschel.esac.esa.int/Docs/SPIRE/html/spire_om.html 
*** For a Gaussian beam the beam power spectrum is also Gaussian, 
and its FWHM F, is related to the half power beam width (HPBW), fl, 

byr= V8/«2/;rD. 



Angular frequency, s [arcmin '] 


Fig.A.l. Power spectrum of the SPIRE 250-/im beam, ybeam)'?) 
(black solid curve). The dashed curve represents the power spec¬ 
trum of a Gaussian beam with a FWHM equal to the nominal 
18'.'2 HPBW of the SPIRE 250 yum beam. The dotted vertical 
line marks the FWHM of the corresponding power spectrum in 
the angular frequency domain. 


umn dens ity map obtained using the multi-scale decomposition 
scheme of iPalmeirim et al.l(l2013ll . The power spectrum values at 
high spatial frequencies has more scatter than the one derived by 
converting I 250 to H 2 -column density (see Fig.|4]l. This is primar¬ 
ily because the multi-scale decomposition technique introduces 
fluctuations at small spatial scales which manifest themselves 
through enhanced scatter in the Fourier modes at high spatial 
frequencies. The distribution of power spectrum slopes derived 
from high-resolution column density maps for the whole sam¬ 
ple of filaments is shown in Fig. 15:21 The mean power spec¬ 
trum slopes before and after correcting for the beam effect are 
ifobs = -1.9 + 0.3 and dcon = -1.6 + 0.45, respectively. Within 
the quoted errors, these results are indistinguishable from those 
found in Sect. 15. B using modified 250 fim maps. 


Angular frequency, s [arcmin '] 

0.1 1.0 4.0 



Spatial frequency, s [pc '] 


Fig.B.l. Power spectrum of the line-mass fluctuations along the 
Pipe filament of Fig. |2]as measured in the high-resolution H 2 - 
column density map resulting from the m ulti-scale decomposi¬ 
tion technique of iPalmeirim et al] (1201 3ll . (Similar to Fig.|?]but 
based on data from the high-resolution column density map in¬ 
stead of the modified 250 fim map.) The power-law fits to the 
power spectra Pobs(s) (black plus symbols) and Ptmeis) (cyan 
filled circles) have logarithmic slopes aobs = -2.3 + 0.2 and 
Qicorr = -1.7 + 0.3, respectively. 
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- 3.0 - 2.5 - 2.0 - 1.5 - 1.0 - 0.5 

Power spectrum slope, a 


Fig.B.2. Distributions of power spectrum slopes measured in the 
high-resolution H 2 -column density maps before beam correction 
(dashed histogram) and after beam correction (solid histogram). 
(Similar to Fig.|7]but based on data from the high-resolution col¬ 
umn density maps instead of the modified 250 jjm maps.) Best-fit 
Gaussian curves to the two observed distributions are overplot¬ 
ted. The two distributions are centered on aobs = -1.9 + 0.3 
and cTcorr = -1.6 + 0.45 for the uncorrected and beam-corrected 
power spectra, respectively. 
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